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The magnetic ground state phase diagram of the ferromagnetic Kondo-lattice model is constructed 
by calculating internal energies of all possible bipartite magnetic configurations of the simple cubic 
lattice explicitly. This is done in one dimension (ID), 2D and 3D for a local moment of S = |. 
By assuming saturation in the local moment system we are able to treat all appearing higher local 
correlation functions within an equation of motion approach exactly. A simple explanation for the 
obtained phase diagram in terms of bandwidth reduction is given. Regions of phase separation are 
determined from the internal energy curves by an explicit Maxwell construction. 



I. INTRODUCTION 

The ferromagnetic Kondo lattice model (FKLM), also 
referred to as s-d model or double exchange model, is 
the basic model for understanding magnetic phenomena 
in systems where local magnetic moments couple ferro- 
magnetically to itinerant carriers. This holds for a wide 
variety of materials. 

In the context of transition metal compounds Zener 
proposed the double exchange mechanism to explain fer- 
romagnetic (FM) metallic phase in the manganites^. In 
these materials the Mn hd shells are split by the crystal 
field into three degenerate t 2g orbitals which are local- 
ized and form a total spin S — | according to atomic 
selection rules and two e g orbitals providing the itinerant 
electrons. These electrons couple via Hund exchange cou- 
pling ferro magnetically with the localized spins. There- 
fore the FKLM is a basic ingredient to describe the rather 
complex physics of the manganites^^. 

Another nearly ideal field of application of the FKLM 
is the description of the rare earth materials Gd and EuX 
(X=0,S,Se,Te). These materials have a half filled 4/ 
shell in common that is strongly localized and the elec- 
trons in this shell couple to a total spin momentum of 
S = | . The FKLM was then used successfully to explain 
the famous redshift of the absorption edge of the opti- 
cal 4:f-5d transition in the ferromagnetic semiconductor 
EuO&£. In | a many-body analysis of the FKLM in 
combination with a band structure calculation was used 
to get a realistic value for the Curie temperature of the 
ferromagnetic metal Gd that is in good agreement with 
experiment. 

Although it is necessary to extent the FKLM in order 
to get a realistic description of the above mentioned ex- 
amples knowledge of the properties of the pure (single 
band) FKLM is crucial for understanding these materi- 
als. 

To reveal the ground state magnetic phases one has 
to solve the many-body problem of the FKLM. This was 
already done in previous works by using different tech- 
niques. Dynamical mean field theory (DMFT) was used 
by several authors [9, ld fli"lfl2| to get information about 
different magnetic domains. In [l3| a continuum field the- 
ory approach was used to derive the 2D phase-diagram 



at T — 0. Classical Monte Carlo simulations were per- 
formed in [9.14]. For ID systems numerical exact density- 
matrix renormalization group calculations were done in 
In [l6[ the authors have used a Green function 
method to test the validity of assuming the quantum lo- 
calized spins to be classical objects. Extended FKLMs 
including more material specific effects were for instance 
investigated in [B1fl7|. 

In this work we will compare all bipartite magnetic 
configurations for the simple cubic (sc) lattice by calcu- 
lating their respective internal energies. To this end the 
electronic Green function has to be determined. This is 
done by an equation of motion approach and, assuming 
that the local moment system is saturated, we are able 
to show that all appearing local higher correlation func- 
tions can be treated exactly. From the calculated internal 
energies the phase-diagram is constructed and region of 
phase-separation are determined. 

The paper is organized as follows. In Sec.|TT]the model 
Hamiltonian and details of the calculation arc presented. 
In Sec. IIIII we discuss the phase-diagrams and give an 
explanation for the sequence of phases obtained by look- 
ing at the quasi-particle density of states. In Sec. IIVI we 
summarize the results and give an outlook on possible 
directions for further research. 



II. MODEL AND THEORY 

A. Model Hamiltonian 

For a proper description of different (anti-) ferromag- 
netic alignments of localized magnetic moments it is use- 
ful to divide the full lattice into two or more sub-lattices 
(primitive cells) each ordering ferro magnetically. 
In this work we only consider simple cubic bipartite lat- 
tices, i.e. anti-ferromagnetic configurations that can be 
obtained by dividing the simple cubic lattice into two 
sub-lattices. In Fig. (JT|) all possible decompositions in 
two and three dimensions are shown. In case of ID only 
the ferromagnetic and g-type anti-ferromagnetic phase 
remain. The Hamiltonian of the FKLM in second quan- 
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FIG. 1: (Color online) Magnetic phases considered in this 
work (ID omitted). 



tization reads as follows: 
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The first term describes the hopping of Bloch electrons 
with spin a between different sites. The lattice sites Ri a 
are denoted by a Latin index i for the unit cell and an 
Greek index a € A, B for the corresponding sub-lattice, 
i.e. Ri Q = R; + r a . The second term describes a 
local Heisenberg-like exchange interaction between the 
itinerant electrons and local magnetic moments Si a 
where J > is the strength of this interaction, z-\[ = ±1 
accounts for the two possible spin projections of the 
electrons and (S[ a = Sf a + z a iS v ia ) denotes the spin 
raising/lowering operator. 



B. internal energy 

The internal energy of the FKLM at T = is given by 
ground state expectation value of the Hamiltonian: 



U=(H) 



f-(E)ES a<7 (E)dE (2) 



where S aa (E) = —^ImG aa (E) is the local spectral den- 
sity, f-(E) denotes the Fermi function and G a<J (E) de- 
notes the local electronic Green function (GF). Note, that 
this formula is obtained by a straightforward calculation 
of the ground-state expectation value of the Hamiltonian 
(P) using the spectral theorem and is therefore exact. 

Our starting point is the equation of motion (EQM) 
for the electronic GF: 
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for the ground state is perfect saturation of the local 
moment system^. With this assumption the Ising-GF 
can be decoupled exactly: 
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where z a = ±1 denotes the direction of sub-lattice mag- 
netization. In a first attempt to solve Eq. ([3]) we have ne- 
glected spin-flip processes completely (F" k "j^ w 0). With 
we then get a closed system of equations which can be 
solved for the electronic GF by Fourier transformation: 
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where e™' 3 is the Fourier transform of the hopping inte- 
gral and a = —a denotes the complementary sub-lattice. 
We will call this solution the "mean-field" (MF) solu- 
tion. Note, that the ferromagnetic phase is contained in 
the above formula by setting e^" to zero. 

To go beyond the MF treatment it is necessary to find 
a better approximation for the spin-flip-GF. To this end 
we write down the EQM for the spin-flip-GF: 
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Our strategy to get an approximate solution for the spin- 
flip-GF is to treat the non-local correlations on a mean- 
field level whereas the local terms will be treated more 
carefully. This is similar to the idea of the dynamical 
mean field theory (DMFT) developed for strongly corre- 
lated electron systems^ Let us start with the non-local 
(i ^ k or i = k but a ^ 7) GFs first. It can be shown-^ 
that the higher GFs resulting from the commutator of 
with H s f are approximately given by the product 
of the spin-flip-GF times spin-wave energies of the local 
moment system. Therefore it is justified to neglect the 
resulting GFs since the spin-wave energies are typically 
3-4 orders of magnitude smaller than the local coupling 
Jim. 

The second term on the rhs of ([6]) gives two higher GFs 
which we decouple on a mean-field level: 
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where in the last step the saturated sub-lattice magneti- 
zation is exploited. 

We now come to the local terms (i = k, a = 7). The two 
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higher GFs resulting from the second commutator on the 
rhs of reduce to: 

((Siu Sia c ia<J ! C j/3 CT )) ~~ * ^(l _ z o-2 Q )G™ CT (8) 

Additionally we get a higher order Ising-GF and spin-flip- 
GF from the first commutator. The higher order spin- 
flip-GF can be treated exactly by using the EQM of the 
(known) Ising-GF given in the appendix (|Al|) . This leads 
to: 

(i^ia n ia<J c iu-o\ c j/3cr)) * 

-{z a z a S-8 aa )F%*. (9) 



The higher order Ising-GF can be traced back to the 
higher order spin-flip-GF by writing down its EQM and 
make use of saturation in the local-moment system (see 
appendix [B] for details): 

-j&,Gs: w «s^-.ii») ■ (io) 

It is a major result of this work that it is possible to 
incorporate all local correlations without approximation, 
i.e. to treat all local higher order GFs exactly. Com- 
bining the results for the appearing higher GFs found 
in 0, ©, © and (JTUJ) we can now solve © for the 
spin-flip-GF: 
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Inserting this result into © and performing a Fourier transformation we finally get: 
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This equation allows for a self-consistent calculation of 
the electronic GF and we will call this the spin-flip (SF) 
solution. 

One important test for the above result is to compare 
it with exact known limiting cases. We found that 
(TT2)) reproduces the solution of the ferro-magnetically 
saturated semiconducto r 21 ! 22 in the limit of zero band- 
occupation. Additionally the 4-peak structure of the 
spectrum as known from the "zero-bandwidth" -limits is 
retained whereas the peaks are broadened to bands with 
their center of gravity at the original peak positions. 




FIG. 2: (Color online) Explicit Maxwell construction for de- 
termining the boundaries of phase separated regions. 



III. RESULTS AND DISCUSSION 



C. phase separation 

To determine the regions of phase separation in the 
phase diagram we have used an explicit Maxwell con- 
struction as shown in Fig[5] The condition for the bound- 
aries of the phase separated region is: 
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The internal energy of the FKLM at T = is given as 
an integral ([2]) over the product of (sub-lattice) quasi- 
particle density of states (QDOS) times energy up to 
Fermi-energy. For understanding the resulting phase- 
diagrams it is therefore useful to have a closer look at the 
QDOS first. In FigEjthe sub-lattice MF-QDOS is shown 
for the different magnetic phases investigated (in 3D). 
The underlying full lattice is of simple cubic type with 
nearest neighbor hopping T chosen such that the band- 
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FIG. 3: Sub-lattice quasi particle density of states (QDOS) of 
up and down electrons obtained from the MF-GF (|6j for two 
values of local coupling J shown for different magnetic con- 
figurations. Parameters: S = | and free electron bandwidth: 
W = 1.0 eV. 



FIG. 4: (Color online) Sub-lattice QDOS of up and down 
electrons obtained from the SF-GF (|12[) for three different 
band-fillings n shown for the ferromagnetic and a-afm phase. 
The local coupling J = 0.5 eV is fixed. Dotted line: corre- 
sponding MF result. Horizontal lines: respective Fermi levels. 
Other parameters as in Fig[3] 



width W is equal to W = 1 eV in the case of free elec- 
trons ( J = eV). The local magnetic moment is equal to 
S = I. We have plotted the up and down-electron spec- 
trum separately for two different values of J = 0.1/1.0 
eV. The exchange splitting A ex — JS eV of up and down- 
band is clearly visible. The decisive difference between 
the phases for nonzero values of J is bandwidth reduc- 
tion from ferromagnetic over a, c to g-afm phase. The 
reason for this behavior becomes clear by looking at the 
magnetic lattices shown in Fig[T] In the ferromagnetic 
case an (up-)electron can move freely in all 3 directions 
of space without paying any additional potential energy. 
In a-type anti-ferromagnetic phase the electron can still 
move freely within a plane but when moving in the di- 
rection perpendicular to the plane it needs to overcome 
an energy-barrier A ex . Hence the QDOS for large values 
of J resembles the form of 2D tight-binding dispersion. 
The bandwidth is reduced due to the confinement of the 
electrons. In the c-afm phase the electron can only move 
freely along one direction and the QDOS becomes effec- 
tively one dimensional. Finally in the g-type phase the 
electron in the large J limit is quasi localized and the 
bandwidth gets very small. We will see soon that this 
bandwidth-effect is mainly responsible for the structure 
of the phase-diagram. Before we come to this point we 
want to discuss the influence of spin-flip processes as in- 
corporated in OH). In FigHthe QDOS for J = 0.5 eV 
is shown for three different band fillings n. The corre- 
sponding Fermi energies are marked by horizontal lines. 
The apparent new feature are the scattering states in 
the down spectrum for band fillings below half filling. 
Thereby the spectral weight of the scattering states is 
more and more reduced with increasing Fermi level. A 
second effect is that the sharp features in the MF-QDOS 
of the anti- ferromagnetic phases are smeared out. Com- 
pared to the MF results the overall change of QDOS be- 
low Fermi energy due to the inclusion of spin-flip pro- 



cesses is small and will not affect the form of the phase- 
diagram drastically. However non-negligible changes can 
be expected. Note that the model shows perfect particle- 
hole symmetry. Therefore the results for the internal en- 
ergy will be the same for n — x and n = 2 — x {x = . . . 1, 
n=l: half filling). 

We come now to the discussion of the phase-diagrams 
which we got by comparing the internal energies of the 
different phases explicitly. The pure phase-diagrams 
(without phase-separation) are shown in Fig|5] whereas 
the different phases are marked by color code. In the 
first column the results of the MF-calculation are shown 
for the 1-, 2- and 3-dimensional case. The second col- 
umn shows the effects of inclusion of spin-flip processes. 
We will concentrate here mainly onto the 3D case since 
most of the given arguments hold equally for the ID 
and 2D case. For J = the system is paramagnetic 
(black bar at bottom). For larger J (J > 0) a typical 
sequence appear: for low band-fillings n the system is 
always ferromagnetic and, with increasing n, it becomes 
a-type then c-type and finally g-type anti-ferromagnetic. 
This behavior is understood easily by looking at the for- 
mula for the internal energy @ and the MF-QDOS in 
FigEl Because of the bandwidth-effect discussed already 
the band-edge of the ferromagnetic state is always low- 
est in energy and will give therefore the lowest inter- 
nal energy for small band-occupation. But since the 
QDOS of the anti-ferromagnetic phases increase much 
more rapidly than the ferromagnetic one, these give more 
weight to low energies in the integral © and will be- 
come lowest in energy eventually for larger band- fillings. 
Therefore the bandwidth-effect is main effect explaining 
the order of phases with increasing n. A very interesting 
feature can be found in the region: J = 0.2 . . . 0.3. In 
this region the ferromagnetic phase is directly followed 
by the c-afm phase for increasing n although the a-afm 
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FIG. 5: (Color online) First column: Phase-diagram as function of band-filling n and local coupling J obtained with mean-field 
(MF) theory © in one, two and three dimensions. Second column: Phase-diagram obtained by inclusion of spin-flip processes 
(SF) (|f 2[) . Regions of different colors mark different (magnetic) phases: ferromagnetic (white), a-afm (brown), c-afm (orange), 
g-afm (green) and paramagnetic phase (black). 



phase has a larger bandwidth than the c-afm phase. This 
can be explained by the two-peak structure of the c-afm- 
QDOS. Due to the first peak at low energies these en- 
ergies are much more weighted than in the a-afm case 
and the c-afm phase will become lower in energy than 
the a-afm phase. Since the reduction of bandwidth of 
the anti-ferromagnetic phases compared to the ferromag- 



netic phase is more pronounced for larger values of J the 
ferromagnetic region is growing in this direction. 
The paramagnetic phase (black bar at J = 0) disappear 
for any finite J since due to the down-shift of the up- 
spectrum of the ferromagnetic phase their internal energy 
will always be lower. When comparing the MF and the 
SF-phase-diagram they appear to be very similar at first 
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FIG. 6: (Color online) Phase-diagram with phase-separation. Regions of phase-separation are marked with two-colored stripes. 
Color code as in Fig[S] 



glance. However two interesting differences can be found, 
namely an increased J region without a-afm-phase and 
the vanishing c-phase above J f» 0.8eV. 

FiglHl shows the phase-diagrams where regions of 
phase-separation, which we have determined by an ex- 
plicit Maxwell construction (|13[) . are marked by colored 
stripes. The two colors denote the involved pure phases. 
As one can see large regions become phase-separated, 
whereas the two participating phases are mostly deter- 



mined by the adjacent pure phases. There is one inter- 
esting exception from this: above a certain J only fm/g- 
afm phase-separation survives and suppresses all other 
phases in this area. Inclusion of spin-flip processes as 
shown in the right column of FigJS] push this J up to 
higher values. Generally spin- flip processes seem to re- 
duce phase-separation as can be seen in the g-afm phase 
and e.g. at the border between fm and c-afm phase. 

Our results are in good qualitative agreement with nu- 
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merical and DMFT results reported by other o 11 ! 12 ' 14 . It 
is common to all these works that for small coupling 
strength J there is only a small ferromagnetic region 
at low band occupation n followed by more complicated 
(anti- ferromagnetic, spiral, canted) spin states/phase- 
separation. With increasing J the region of fm is also 
increased to larger n values. Near half-filling (n = 1) one 
will find always anti-ferromagnetism/phase-separation. 
phase-diagram very similar to our 2D-FM result shown in 
FigJBjwas obtained by Pekker et.al.JHj]. The positions of 
A and G phase are in nearly perfect agreement. However 
the authors seem not to have taken into account phase- 
separation between A and G phase and their finding of 
FM/A phase-separation near half- filling at larger J is not 
in accordance with our results. 



when a saturated sub-lattice magnetization is assumed. 
The EQM of the remaining GF turns out to be: 
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IV. SUMMARY AND OUTLOOK 



Subtracting the term denoted by (I) from this equation 
one gets: 



We have constructed phase diagrams of the FKLM in 
ID, 2D and 3D by comparing the internal energies of all 
possible bipartite magnetic configurations of the simple 
cubic lattice. To this end the electronic GF is calculated 
by an EQM approach. We can show, that it is possible 
to treat all appearing higher local correlation functions 
exact and we derive an explicit formula for the electronic 
GF (fT2")l . The obtained sequence of phases with increas- 
ing band occupation n and Hunds coupling J is explained 
by the reduction of QDOS bandwidth due to electron 
confinement. Region of phase separation are then de- 
termined from the internal energy curves by an explicit 
Maxwell construction. 

In the phase diagram obtained only phases appear that 
have explicitly considered by us. Therefore an important 
extension of this work could be the inclusion of more 
complicated spin structures like canted/spiral spin states 



as reported by others 13lfl5l |. However the bandwidth 



criterion obtained here can certainly be applied to such 
more complicated states also. 
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This can be solved for ((ri,j a _ CT Ci acr ; c^ CT )) by left- 
multiplying with the MF-GF matrix: 



APPENDIX A: EQM OF THE ISING-GF 
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APPENDIX B: HIGHER ORDER ISING-GF 
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The higher order Ising-GF can be decomposed into: 



Two other equations are obtained from (|B2p by subtract 



((S^riia-aCtaa-; c Jj3(J )) — > z a S {{niot-vCitxo) c j0 a )) C^l) m S term (II) or (III) and performing the same steps as 
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before. This yields: Adding (|B5|) and (IB6[) and subtracting (|B4[) one finally 

gets: 
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